Fluctuations in network dynamics 
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Most complex networks serve as conduits for various dynamical processes, ranging from mass 
transfer by chemical reactions in the cell to packet transfer on the Internet. We collected data 
on the time dependent activity of five natural and technological networks, finding that for each 
the coupling of the flux fluctuations with the total flux on individual nodes obeys a unique scaling 
law. We show that the observed scaling can explain the competition between the system's internal 
collective dynamics and changes in the external environment, allowing us to predict the relevant 
scaling exponents. 
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Recent advances in uncovering the mechanisms shap- 
ing the topology of complex networks are overshad- 
owed by our lack of understanding of common organizing 
principles governing network dynamics. In particular, we 
are far from understanding how the collective behavior 
of often millions of nodes contribute to the observable 
dynamical features of a given system, prompting us to 
continue the search for dynamical organizing principles 
that are common to a wide range of complex systems. 
To make advances in this direction we need to comple- 
ment the available network maps with data on the time 
resolved activity of each node and link. 

Traditional approaches to complex dynamical systems 
focus on the long time behavior of at most a few dynam- 
ical variables, characterizing either a single node or the 
system's average behavior. To simultaneously character- 
ize the dynamics of thousands of nodes we investigate the 
coupling between the average flux and fluctuations. Our 
measurements indicate that in complex networks there 
is a characteristic coupling between the average flux (/*) 
and dispersion m of individual nodes (Fig. 1). To quan- 
tify this observation we plot <Tj for each node i in func- 
tion of the average flux (/.;} of the same node (Figs. 
|3J). We find that for five systems for which extensive dy- 
namical data is available the dispersion depends on the 
average flux as 
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Most intriguing, however, is the finding that the dynam- 
ical exponent a is in the vicinity of two distinct values, 
a = 1/2 (Fig. and a = 1 (Fig. EJ, suggesting that 
diverse real systems can display two distinct dynamical 
universality classes. 

The a ~ 1/2 systems (Fig. 2): The Internet, viewed 
as a network of routers linked by physical connections, 
serves as a transportation network for information, car- 
ried in form of packets Q • Daily traffic measurements 
of 374 geographically distinct routers indicate that the 
relationship between traffic and dispersion follows for 
close to seven orders of magnitude with a 1 = 1/2 (Fig. 

In a microprocessor, in which the connections be- 
tween logic gates generate a static network Q , informa- 
tion is carried by electric currents. At each clock cycle 
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FIG. 1: (a) Time dependent traffic on three Internet routers 
of the Mid- Atlantic Crossroads network, whose activity is 
monitored by the Multi Router Traffic Grapher software 
(MRTG). The figure shows the number of bytes per second 
for each of the routers in five minutes intervals for a two day 
period, (b) Streamflow, measured in cubic feet per second, 
on three rivers in the US river basin, based on data collected 
by the U.S. Geological Survey in 2001. On the right of each 
plot we show the time average of the flux (/) displayed as 
horizontal dotted lines superposed on the graphs, and the 
dispersion, a, for each signal, indicating orders of magnitude 
differences in both flux and dispersion between nodes of the 
same network. 



a certain subset of connections i are active, the relevant 
dynamical variable fi{t) taking two possible values, or 
1. The activity during 8,862 clock cycles on 462 nodes 
of the Simplel2 microprocessor indicates that the aver- 
age flux and fluctuations follow (J), with a m — 1/2 (Fig. 

03). 

The a ~ 1 systems (Fig. 3): The WWW, an exten- 
sive information depository, is a network of documents 
linked by URLs [4|. As many websites record individ- 
ual visits, surfers collectively contribute to a dynamical 
variable f™(t) that represents the number of visits site i 
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FIG. 2: The relationship between fluctuations (a) and the av- 
erage flux ({/}) for the a = 1/2 systems, (a) Time resolved 
information for 374 Internet routers of the Mid- Atlantic 
Crossroads, ABILENE network, MIT routers, UNAM routers, 
all Brazilian RNP backbones, and dozens of smaller routers 
on the Internet, covering for each node two days of activity 
with five minute resolution, (b) The activity of the 462 signal 
carriers of the 12-bit Simplel2 microprocessor, recorded over 
8, 862 clock cycles. 
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FIG. 3: The relationship between fluctuations (<r) and the 
average flux ((/)) for systems belonging to the a = 1 class, 
(a) Daily visitations on websites collected using the Nedstat 
web monitor. We analyzed daily traffic for a 30 day period 
for 1,000 sites in USA (circles) Brazil (squares) and Japan 
(triangles), (b) The daily streamflow of 3,945 rivers on the 
US river basin during the year of 2001 is recorded by the US 
Geological Survey, (c): Daily traffic on Colorado and Ver- 
mont highways representing the daily number of cars pass- 
ing through observation points on 127 highways from 1998 to 
2001. 



receives during day t. We studied the daily breakdown of 
visitation for 30 days for 3, 000 sites scattered over three 
continents, determining for each node i the average (ff) 
and dispersion of. As Fig. shows, of and (ff) fol- 
low 0J over three orders of magnitude with dynamical 
exponent a w = 1. The highway system is an example of 
a transportation network, the relevant dynamical vari- 
able being the traffic at different locations. We analyzed 
the daily breakdown of traffic measurements at 127 loca- 
tions on Colorado and Vermont highways. The results, 
shown in Fig. again document scaling spanning over 
five orders of magnitude with a — 1. Finally, the river 
network is a natural transportation system R whose dy- 
namics is probed via time resolved measurements on the 
stream of several US rivers on 3,495 different locations. 
While these fluctuations are driven by weather patterns, 
the relationship between the average stream and its fluc- 
tuations again follows with a r = 1 (Fig. |2t). 

To understand the origin of the observed dynamical 
scaling law Q we study a simple dynamical model that 
incorporates some key elements of the studied systems. 
While the topology of these systems vary widely, from a 
tree (rivers) to a scale-free network (WWW, Internet), a 
common feature of the studied systems is the existence of 
a transportation network that channels the flux toward 
selected nodes. Therefore, we start with a network of 
N nodes and L links, described by an adjacency matrix 
Mij, which we choose to describe either a scale-free or 
a random network 0. As the dynamics of the studied 
systems varies widely, we study two different dynamical 
rules. Model 1 considers the random diffusion of W walk- 
ers on the network, such that each walker that reaches 
a node i departs in the next time step along one of the 
links the node has. Originally each walker is placed on 
the network at a randomly chosen location and removed 
after it performs M steps, mimicking in a highly sim- 
plified fashion a human browser surfing the Web for in- 



formation. To probe the collective transport dynamics 
counters attached to each node record the number of vis- 
its by various walkers. To capture the day to day fluc- 
tuations on individual nodes we repeat independently D 
times the diffusion of W walkers on the same fixed net- 
work and denote by fi(t) the number of visits to node i 
on day t = 1, . . . , D. As Fig. 0Ji indicates, the average 
flux and fluctuations follow (JJJ with a = 1/2. In Model 
2 we replaced the diffusive dynamics with a directed flow 
process. In this case each day t we pick W randomly 
selected pairs of nodes, designating one node as a sender 
and the other as a recipient, and send a message between 
them along the shortest path. Counters placed on ev- 
ery node count the number of messages passing through. 
This dynamics mimics, in a highly schematic fashion, the 
low density traffic between two nodes on the Internet. 
As Fig. 4d shows, we find that Model 2 also predicts 
q = 1/2, indicating that the a — 1/2 exponent is not a 
particular property of the random diffusion model, but it 
is shared by several dynamical rules. 

We can understand the origin of the a — 1/2 expo- 
nent if we inspect the nature of fluctuations in Model 1. 
In the M = 1 limit walkers arrive to randomly selected 
nodes but fail to diffuse further, reducing the dynam- 
ics to random deposition, a well known model of surface 
roughening [6( . Therefore, the average visitation on each 
node grows linearly with time, (/) ~ t, and the disper- 
sion increases as o ~ i 1/2 , providing a = 1/2 @. While 
for M > 1 diffusion generates correlations between the 
nodes, we find that the fluctuations on the individual 
nodes, <t™*, continue to be dominated by the internal 
randomness of the walker arrival and diffusion process, 
following the a = 1/2 dynamical exponent [7j . 

To understand the origin of the second (a = 1) univer- 
sality class, we note that in real systems the fluctuations 
on a given node are determined not only by the system's 
internal dynamics, but also by changes in the external 
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environment. To incorporate externally induced fluctua- 
tions we allow W (the number of walkers and messages 
in Models 1 and 2), to vary from one day to the other. 
Assuming that the day to day variations of W(t) define a 
dynamic variable chosen from an uniform distribution in 
the interval [W- AW, W + AW], for AW = we recover 
a = 1/2. However, when AW exceeds a certain thresh- 
old, in both models the dynamical exponent changes to 
a = 1 (Fig. 2t> and e). 




FIG. 4: In Model 1 on each "day" t we release W(t) = 
(W) + £(t) walkers on randomly selected nodes and allow 
them to perform M = 10 3 random diffusive steps, where £(t) 
is a uniformly distributed random variable between — AW 
and AW and (W) = 10 4 . (a) The figure shows the a ((f)) 
curves for AW = 0, 20, 40, 80, 100, 200, 800, 1000, 4000, 10000 
from top to bottom, (b) The dependence of the exponent a 
on AW, obtained by fitting the <n versus (fi) curves shown 
in (a). Note that while the figure shows a gradual transi- 
tion, the transition in infinite systems should be sharp, (c) 
Average fluctuations (<Ti), obtained by averaging Oi over all 
nodes i in the system, shown in function of the amplitude 
of the external driving force AW. While under AW ~ 10 2 
the magnitude of (<7j) is independent of AW, for large AW 
the fluctuations increase rapidly, indicating that the network 
dynamics is externally driven, (d-f) The same as in (a-c), 
but for Model 2, where the diffusive dynamics was replaced 
by message passing. W was again chosen from an uniform 
distribution of width AW and average (W) = 10 4 . In all sim- 
ulations we used a scale-free network [3 with 7 = 3 and 10 4 
nodes. 

To understand the origin of the a = 1 exponent we 
notice that on each node the observed day to day fluc- 
tuations have two sources. For AW = we have only 
internal fluctuations, coming from the fact that under 
random diffusion (or random selection of senders and re- 
ceivers in Model 2) the number of walkers (messages) 
that visit a certain node displays day to day fluctua- 
tions. For AW ^ the fluctuations have an external 
component as well, as when the total number of walkers 
(messages) change from one day to the other, they pro- 



portionally alter the visitation of the individual nodes 
as well. If the magnitude of the day to day fluctuations 
is significant, they can overshadow the internal fluctu- 
ations cr] n *. Indeed, if in a given time frame the to- 
tal number of walkers or messages doubles, the flux on 
each node is expected to grow proportionally, a poten- 
tially much larger variation than the changes induced by 
the internal fluctuations. Therefore, for AW 7^ the 
external driving force, determined by the time depen- 
dent W(t), contributes to the daily fluctuations with a 

dispersion a dr (AW) = ^ (W(t) 2 ) - (W(t)) 2 . The to- 
tal fluctuations for node i are therefore given by of = 

(o~l nt ) + (af xt ) 2 . As the effect of the driving force is 
felt to a different degree on each node, we can write 
a ext _ Aia dr (AW), where Ai is a geometric factor cap- 
turing the fraction of walkers channeled to node i, and 
depends only on the position of node i within the net- 
work. When AW = 0, the external component a dr van- 
ishes, resulting in er™* = cii (fi) 1 ^ 2 , as discussed earlier, 
where a, is an empirically determined coefficient. When 
AW is sufficiently large, so that A 4 cr dr (AIT) > erf 1 *, 
then the fluctuations on each node are dominated by 
the changes in the external driving force. In this limit 
a node's dynamical activity mimics the changes in the 
external driving force, allowing us to approximate the 
flux at node i with /j(t) = AiW(t). In this case we 

have (fi) = A, (W(t)) and (f 2 ) = A 2 (w(t) 2 ), giving 

°i = yj (U 2 ) - {fif = A t a dr . As a dr and (W(t)) are 
time independent characteristics of the external driving 



force, we find <7j ~ af xt — 
served coupling with a 



(fi), providing the ob- 
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= 1. Note that this derivation 
is independent of the network topology or the transport 
process, predicting that any system for which the magni- 
tude of fluctuations in the external driving force exceeds 
the internal fluctuations will be characterized by an a = 1 
exponent. 

These calculations imply that the fluctuations on a 
given node can be decomposed into an internal and an 
external component as 



(fi) 



(fi) 
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Therefore, increasing the amplitude of fluctuations AW 
should induce a change from the a = 1/2 intrinsic or en- 
dogenous to the a = 1 driven behavior. To confirm the 
validity of this prediction, in Figs. 4c and f we show the 
average fluctuation over all nodes in function of the 
amplitude AW of the driving force. For both models we 
find that for small AW values di remains unchanged, as 
in this regime di ~ a} nt > <rf xt . However, after AW ex- 
ceeds a certain threshold, di changes behavior, monotoni- 
cally increasing with AW. In this second regime the fluc- 
tuations are driven by external forces, di ~ crf xt ~ Ai<r dr ', 
and according to J2J) we should observe a = 1 . Indeed, we 
find that in both models the transition from the constant 
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to the increasing <fj regime (Figs. 4c, f) coincides with the 
crossover from the a = 1/2 to a = 1 (Figs. 4b, e). Note, 
however, that the gradual transition observed in Figs. 
4b-e from a = 1/2 to a — 1 is a numerical artifact of the 
fitting process: in the transition regime the a — 1/2 and 
a = 1 scaling coexist on the same cr((/)) curve, giving an 
exponent that is different from 1/2 or 1. In reality the 
transition between the two regimes is sharp. To under- 
stand to what degree our findings depend on the specific 
simulation and model details we changed the topology 
from scale-free to random network and from undi- 
rected to directed network, as well as altering the nature 
of the external fluctuations by keeping W constant in 
Model 1 but forcing the number of steps, M, to play the 
role of the stochastic external driving force. For each ver- 
sion we recover the transition between the a = 1/2 and 
a = 1 when the amplitude of the external fluctuations 
exceeds a certain threshold 0] • 

These results indicate that the a — 1/2 exponent cap- 
tures an endogenous behavior, determined by the sys- 
tem's internal collective fluctuations. In the studied 
model internal fluctuations are rooted in the randomness 
in the walkers' arrival and diffusion; on the Internet they 
originate in the choices users make to where and when to 
send a message; for the computer chip they come from the 
alternating utilization of the various circuits, as required 



by the performed computation. In contrast, the a = 1 
exponent describes driven systems, in which the fluctua- 
tions of individual nodes are dominated by time depen- 
dent changes in the external driving forces. Therefore, 
fluctuations of World Wide Web traffic, river streams 
and highway traffic are driven by such external factors 
as daily variations in the number of Web surfers, sea- 
sonal or daily changes in precipitation or daily variations 
in the number of drivers, respectively. 

Of the two observed exponents our derivation indicates 
that a = 1 is universal, being independent of the na- 
ture of the internal dynamics or the network topology. 
There are no firm restrictions, however, on the scaling 
of the internal dynamics, raising the possibility that self- 
organized processes could lead to collective fluctuations 
that are characterized by a exponents different from 1/2. 
Empirical evidence for potential intermediate a values 
comes from ecology, where Q describes spatial and tem- 
poral variations of populations 0] . It is much debated, 
however, whether the observed scaling represent valid ex- 
ponents, or only crossovers between a = 1/2 and 1 |12| . 

We are indebted to Jay Brockman and Steven Balen- 
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research was supported by grants from NSF, NIH and 
DOE. 



[1] S. Bornholdt, H.G. Schuster, Eds. Handbook of Graphs 
and Networks (Wiley- VCH, Berlin, 2002); S.N. Dorogovt- 
sev, J.F.F. Mendes, Evolution of Networks: From Biolog- 
ical Nets to the Internet and WWW (Oxford University 
Press, Oxford, 2003); R. Albert, A.-L. Barabasi, Rev. 
Mod. Phys. 74, 47 (2002); S.H. Strogatz, Nature 410 268 
(2001). 

[2] A. Vazquez, R. Pastor-Satorras, A. Vespignani, Phys. 

Rev. E 65, 066130 (2002). 
[3] R. F. Cancho, C. Janssen, R.V. Sole, Phys. Rev E, 64, 

046119 (2001). 
[4] S. Lawrence, L. Giles, Science 280, 98 (1998). 
[5] J.R. Banavar, A. Maritan, A. Rinaldo, Nature 399, 130- 

132 (1999); M. Cieplak, A. Giacometti, A. Maritan, A. 

Rinaldo, J.R. Banavar, J. Stat. Phys. 91, 1 (1998); G. 

Caldarelli, Phys. Rev. E 63, 21118 (2001). 
[6] A.-L. Barabasi, H.E. Stanley, Fractal Concepts in Surface 

Growth (Cambridge University Press, Cambridge, 1995); 

F. Family, T. Vicsek, Eds., Dynamics of Fractal Surfaces 

(World Scientific, Singapore, 1991). 
[7] While Poissonian statistics explains the a = 1/2 scal- 



ing in the random deposition model, it is not appropri- 
ate for Ethernet traffic |g|, indicating that the Internet's 
a — 1/2 exponent is rooted in the interplay between net- 
work topology and dynamics. The exponential distribu- 
tion describing waiting times is known to predict a = 1. 
[8] W.E. Leland, M.S. Taqqu, W. Willinger and D.V. Wil- 
son, IEEE/ACM Trans. Net., 2, No. 1, 1 (1994). 
[9] A.-L. Barabasi and R. Albert, Science 286, 509 (1999). 

[10] The relative strength of fluctuations can be determined 
independently of the scaling exponents. Such measure- 
ments indicate that for the Internet and the microchip 
the magnitude of the internal fluctuations is close to 100 
times larger than the magnitude of the external fluctua- 
tions. In contrast, for the other systems the internal and 
external fluctuations are comparable in magnitude [M. 
Argollo de Menezes and A.-L. Barabasi, to be published]. 

[11] L.R. Taylor, Nature 189, 732 (1961). 

[12] R.M. Anderson, D.M. Gordon, M.J. Crawley, M.P. Has- 
sel, Nature 296, 245 (1982). 



